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Abstract 

Many of the interesting patterns seen in recent multi- frequency Faraday experiments 
can be understood on the basis of three- wave interactions (resonant triads). In this 
paper we consider two-frequency forcing and focus on a resonant triad that occurs 
near the bicritical point where two pattern-forming modes with distinct wavenum- 
bers emerge simultaneously. This triad has been observed directly (in the form of 
rhomboids) and has also been implicated in the formation of quasipatterns and su- 
perlattices. We show how the symmetries of the undamped unforced problem (time 
translation, time reversal, and Hamiltonian structure) can be used, when the damp- 
ing is weak, to obtain general scaling laws and additional qualitative properties of 
the normal form coefficients governing the pattern selection process near onset; such 
features help to explain why this particular triad is seen only for certain "low" forc- 
ing ratios, and predict the existence of drifting solutions and heteroclinic cycles. We 
confirm the anticipated parameter dependence of the coefficients and investigate its 
dynamical consequences using coefficients derived numerically from a quasipotential 
formulation of the Faraday problem due to Zhang and Vihals. 
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1 Introduction 



When a container of fluid is shaken vertically with sufficient strength, patterns of stand- 
ing waves (SW) develop on the free surface. This well-known transition, first described by 
Faraday [1], is one of the best known examples of a pattern forming instability. At a criti- 
cal value of the forcing amplitude the flat surface loses stability and, with the aid of small 
perturbations, gives way to a spatially structured state. This excited state is often highly 
ordered and has been observed in the form of rolls, squares, hexagons, and more exotic struc- 
tures such as superlattices and quasi-patterns having 8, 10, or 12-fold rotational symmetry 
[2,3,4,5,6]. The temporal spectrum of these resonant SW oscillations is highly ordered as 
well. In the standard Faraday experiment with single frequency forcing of period T = 2tt/uj 
the first pattern to appear is a subharmonic response, i.e., it is 2T-periodic in time [7]. The 
characteristic wavenumber k of the pattern is determined by the dispersion relation Q(k) 
for gravity-capillary waves through the subharmonic resonance condition Q(k) = uo/2. This 
frequency u, which controls the spatial length scale, and the forcing amplitude / are the 
most readily accessible control parameters. 

Recently, a number of investigators [4,5,6,8,9,10,11,12] have looked at more general kinds 
of periodic forcing. Such forcing functions have multiple frequency components and present 
a larger number of adjustable parameters. One of the biggest advantages to using multiple 
frequencies is that one can study, in a controlled manner, the competition between modes 
of different length scales. If the applied acceleration f(t) is composed of two commensurate 
frequencies of ratio m:n (m and n coprime): 

/(*) = \ (fme mu;t + f n e mut + c.c.) = |/ m | cos(mut + <j> m ) + |/ n | cos(n^ + n ), (1) 

an appropriate balance between |/ m | and |/ n | ensures that the two modes (one driven primar- 
ily by the f m forcing component and the other by /„) onset simultaneously. This bicritical 
point is the natural organizing center for the two-frequency problem because it is at this 
point that the greatest variety of patterns evolve and compete with each other while still 
in the weakly nonlinear regime. A complete description of this competition is exceedingly 
difficult to achieve, yet a number of simple resonant interactions have emerged as crucial 
building blocks for many of the more interesting patterns. In this respect resonant triads 
are particularly influential because they are the lowest order nonlinear effect (i.e., they arise 
at second order in a weakly nonlinear expansion) and provide a strong phase coupling be- 
tween the three waves involved. They have been implicated in a number of pattern formation 
studies, both experimental [3,4,5,8,13] and theoretical [10,14,15,16,17]. 

In this paper we consider the two-frequency forcing of Eq. (1) and focus on a simple three- 
wave interaction, occurring near the bicritical point, in which two wavevectors from one 
critical circle, k\ and &2, sum to give a third wavevector on the second critical circle 
(see Fig. 1). Patterns formed by these three SW modes were called two-mode rhomboids by 
Arbell and Fineberg [8]. The same basic resonance, however, was shown to be an essential 
component of more complicated states such as quasipatterns with 8-fold or 10-fold symme- 
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Fig. 1. Resonant triad interactions at the bicritical point: (a) when < k\ and (b) 
when k\ < k% < 2k±. 

try [8] (the former observed with three- frequency forcing [3]), and several other multiple 
mode and modulated states [3]. One of the most striking facts about this resonant triad is 
that it has so far been seen only with forcing ratios of 1:2 (as part of a superlattice-II type 
state [3]), 2:3, and 4:5. In [9] we showed that the absence of two-mode rhomboids for larger 
values of m and n is related to the presence of (broken) time translation symmetry. This 
paper provides additional details behind that result as well as others in [9], and explores 
their dynamical implications. 

Symmetry arguments, in general, play a central role in our understanding of pattern for- 
mation, providing a powerful framework that explains generic features, i.e., properties that 
one expects to find across a wide range of physical systems sharing the same symmetry. 
For Faraday systems of large aspect ratio (i.e., where the horizontal extent of the system 
is large compared to the characteristic wavelength) the relevant symmetry is that of the 
plane, described by the Euclidean group E(2). If the pattern is spatially periodic, as is often 
the case near onset, many useful results (see, e.g., [18,19]) can be obtained by restricting to 
patterns commensurate with this periodic lattice. The choice of a lattice fixes the orientation 
and removes the arbitrary rotations that are part of E(2). Spatial translations, however, are 
still allowed, as are certain discrete rotations (60° for hexagons, for example) and reflections. 
Thus, when periodic planforms are considered the noncompact group E(2) is replaced by a 
compact group, namely, the semi-direct product of a discrete group representing the lattice 
symmetries (this group is D 2 for the resonant triads we consider) and the two-torus rep- 
resenting spatial translations. The fact that these symmetry operations must transform a 
given solution into an equivalent solution restricts the form of the equations governing the 
pattern formation process near onset and thus determines which types of nonlinear terms 
can be present and which cannot [18]. 

Temporal symmetries also play an important role in periodically forced systems. In par- 
ticular, such systems are invariant under discrete translation through one forcing period: 
t — > t + T. Whether this symmetry is important or not depends on the nature of the linear 
eigenfunctions. In the case of single-frequency forcing, for example, the eigenfunctions are 
subharmonic (i.e., they have odd parity under t — > t + T) and the associated Z2 symmetry 
bars all even terms from the evolution equations; in particular, the absence of quadratic 
terms means that there are no resonant triad interactions. With two-frequency forcing the 
restrictions due to discrete temporal translation are less severe because the period T is longer 
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than the corresponding period for, say, isolated moo forcing, which would be T/m. In this 
way, modes driven "subharmonically" by f m may actually be harmonic with respect to the 
full period T; they are harmonic if m is even and subharmonic if m is odd, similarly for n. 

In this paper we follow the approach outlined in [9] and extend the usual analysis based on 
"exact" spatial symmetries and discrete time translation to include the approximate (i.e., 
weakly broken) symmetries of continuous time translation, time reversal, and Hamiltonian 
structure. In [9] we considered hexagons and two-mode superlattices [3,5] as well as two- 
mode rhomboids, and showed how these symmetries can be used to determine the leading 
order dependence of important coefficients in the SW amplitude equations. In this paper 
we supply the details of that calculation for the case of two-mode rhomboids (here simply 
referred to as resonant triads) and explore the dynamical consequences of its predictions. 
For example, in Section 4 we show that a preference for resonant coefficients of opposite 
sign [9], which stems from the presence of Hamiltonian structure in the undamped problem, 
leads to the existence of drifting solutions, modulated waves, and heteroclinic cycles. 

The organization of the paper is as follows. In Section 2 we demonstrate how the broken 
temporal symmetries can be replaced by unbroken parameter symmetries (see also [9]), and 
use these in conjunction with the spatial symmetries to determine not only the form of the 
SW amplitude equations but also the manner in which the coefficients of these equations 
depend on the forcing parameters (m, n, |/ m |, |/ n |, <f) m , <f) n ) and on the damping. In Section 3 
we test our symmetry-based predictions against coefficients calculated numerically from a 
quasipotential formulation of the Faraday problem due to Zhang and Vihals [14,15]. Section 4 
summarizes the typical dynamics associated with the resonant triad equations and uses 
numerical examples to illustrate the effects of the broken symmetries discussed in Section 2. 
In Section 5 we offer our conclusions. 



2 Derivation of the Amplitude Equations 

In this section we use general symmetry arguments to derive the form of the equations 
governing resonant triads in the presence of two- frequency forcing (Eq. (1)). In particular, 
we obtain general scaling properties of the normal form coefficients. These scaling laws rest 
on the assumption of weak dissipation (the critical forcing amplitude needed to destabilize 
the spatially uniform state is then similarly small). In this limit the problem is constrained 
not only by spatial symmetries but by the (weakly broken) symmetries of time translation 
t — > £ + s, sGM and time reversal t — > — t that are features of the unforced inviscid 
hydrodynamical problem. In fact, if these two symmetries are extended in an appropriate way 
(see Section 2.2) they apply to the damped forced problem as well. Additional constraints 
resulting from full Hamiltonian structure in the undamped forced problem are examined in 
Section 2.3.1. 

There are many adjustable parameters in the two-frequency Faraday experiment. In addition 
to the fluid depth (assumed to be infinite in our numerical study; see Section 3) and proper- 
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ties such as viscosity v, density p, and surface tension T, one can vary the basic frequency uj 
of the driving, the ratio of the two frequency components m:n, their respective magnitudes 
\f m \ and \f n \, and their phases <p m and <p n . Here we take 1/uj as the unit of time, writing 
r = cut, and use the wavenumber k associated with Q = min(m, n)cu/2 via the inviscid 
dispersion relation goko + Tk^/ p = Q% (g is the acceleration of gravity) as the inverse unit of 
length. This scaling groups cu, v, T, and p into three nondimensional parameters: 7 = 2vk\juj 
(the damping), To = Tk^/lpuj 2 ) (the capillarity number), and Go = goko/u 2 (the gravity 
number). The last two are related by the equation r + G = (min(m, n)/2) 2 . 

The situation of special interest here is the bicritical point where two modes of distinct 
wavenumber are simultaneously excited as the overall forcing amplitude / = (|/ m | 2 + \fn\ 2 ) 1 ^ 2 
is increased (see Fig. 2). We think of this bistability criteria as determining the ratio |/ m |/|/n|- 
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Fig. 2. Neutral stability curves at the bicritical point calculated from the 
Zhang- Virials equations (see Section 3). We use m:n = 3:2, 7 = 0.1, Tq = 0.5, 
G = 0.5, — 0n — and I/2I/I/3I — 0.423. The applied acceleration / is in units 
of go- 

The two critical circles at the bicritical point and the two different ways of constructing a 
resonant triad are illustrated in Fig. 1. In either case we choose three wavevectors: two of 
equal magnitude, |fei| = \k%\ = k\, and their sum, &3 = fci + &2. Without loss of generality 
we associate k± and k 2 with the f m forcing component and fc 3 with the /„ component. Thus 
Fig. la applies if m > n while Fig. lb applies if m < n. As mentioned, when the damping is 
small each of the excited modes can, to first approximation, be thought of as a subharmonic 
response to the corresponding component (f m or f n ) of the driving. The expjifci • x} and 
exp{i&2 • x} modes oscillate with a dominant frequency of m/2 while the exp{ik^ ■ x} mode 
has dominant frequency n/2. A given mode is harmonic (subharmonic) when m or n is even 
(odd). As shown in Section 2.2, resonant triads only occur if n is even, in which case m is 
odd since m and n are coprime. 

As long as damping is present the first modes to appear as / is increased are SW. Sufficiently 
close to onset the resonant triad system may therefore be described by the three complex 
amplitudes Aj, j — 1, 2, 3 of the SW modes exp{ikj ■ x}. We do not begin immediately with 
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these three SW modes but, following the program of [9], obtain first the form of the equations 
governing six traveling wave (TW) modes having the same wavevectors k\, k 2 , and fe 3 . It 
should be noted that in the absence of dissipation and forcing the Faraday system supports 
TW at infinitely many frequencies and wavenumbers, i.e., the only requirement is from the 
inviscid dispersion relation Q(k). It is precisely the applied forcing which selects, through the 
temporal resonance condition, the frequencies and hence the length scales of the problem. 
Our approach, by restricting attention to the wavenumbers k\ and k^, anticipates the effect 
of the parametric forcing. At the same time, however, it assumes that the amplitude of 
this forcing (and of the damping) is small. In this limit a TW expansion is justified despite 
the fact that pure TW solutions are destroyed by the parametric forcing term - such an 
approach may also be viewed as the unfolding of a Hopf bifurcation (occurring at 7 = 0) 
with nearly resonant temporal forcing (see [20,21]). Here we are motivated in our use of TW 
by the action of the weakly broken temporal symmetries; TW transform in a simple way 
under time translation r — > r + s and time reversal r — > — r whereas SW do not. The simple 
transformation properties of TW allow us to take full advantage of the spatial and temporal 
symmetries before a center manifold reduction is used to obtain the equations describing 
the three SW amplitudes A y 

We now systematically consider the two unfolding parameters which move the system away 
from the "ideal" undamped, unforced case: the damping parameter 7 > (proportional to 
viscosity) that breaks time reversal symmetry, and the forcing function /(r) (equivalently, 
f m and f n ), which breaks the continuous time translation symmetry and, in general, the 
time reversal symmetry as well. In the process we pay special attention to resonant terms, 
i.e., those terms which couple the phases of the different modes, because of the crucial role 
they play in the ensuing dynamics. The leading resonant terms are quadratic in amplitude 
but, with the exception of m:n = 1:2, appear only as a result of broken time translation 
symmetry. The residual effect of this symmetry determines, for small 7, the influence (i.e., 
size) of these resonant terms by prescribing how they scale with the symmetry-breaking 
parameters f m and f n . 



2.1 Traveling wave equations 



We consider a system of six undamped, unforced TW with wavevectors fci, fc 2 , and fc 3 . For 
example, the surface height h(x, r) is written as 

-J-) — ^+g»(fcl-!B+mT/2) _|_ ^- e i{k\-x-mT /2) _|_ ^+ ^i(k 2 -x+mT /2) _|_ r^- ^,{k 2 -x-rxiT j2) 

z + e i(k 3 . X +nr/2) + z - e i{k 3 - X -nr/2) + Q C + _ ^ 

The equations describing the slow evolution of the TW amplitudes Zf, j = 1,2,3 must 
be equivariant under the operations of spatial translation (%/,), reflection (k) about a line 
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through k 3 , and rotation (1Z) by n (i.e., inversion through the origin): 



T : Zf -> Z±e^, = (0 b 2 ), 0!, 2 G [0, 2tt), 3 = 0i + 2 , (3) 
k : Z± <- Z±, (4) 
• Zf. (5) 

In the absence of forcing (/) the equations must also be equivariant under time translation: 
T s : (Z± Z±) -> (Z± Z 2 ± )e ± ^/ 2 , Z± - Zfe ±ms '\ s G [0, 4tt), (6) 

while if both / and 7 are zero there is the time reversal symmetry a: 

a :r^-r, Z f - Zf , (7) 

The constraints of symmetries (3-7) lead to equations, truncated at cubic order, of the form 

Z+ = {aZ+Z+} 1:2 + bZ+Z^Z^ 

+ (a|Z+| 2 + b|Z~| 2 + c|Z+| 2 + d|Z~| 2 + e|Z+| 2 + f |Z 3 ~| 2 ) Z+, (8) 



Zt = {/3Z+Z+} 1:2 + (g(|Z+| 2 + |Z+| 2 ) + h(|Zr| 2 + |Z 2 1 2 ) + 1|Z 3 + | 2 + p|Z 3 -| 2 ) Z 3 + , (9) 

where the overdot indicates a derivative with respect to r. The other four equations follow 
from the symmetries k and 1Z. The quadratic terms in curly brackets (subscript 1:2) are 
permitted only when m:n = 1:2 because with other forcing ratios they violate the time 
translation symmetry T s . Whether these resonant terms are present or not, the time reversal 
symmetry a (more precisely, a o 1Z) forces all coefficients to be purely imaginary: 

a = i&i, 5 = i5i, a = ia>i, p = ipi. (10) 

The absence of linear terms with imaginary coefficients is not related to symmetry but to 
the way we define k\ and k 3 , which requires that the detuning vanish at / = 7 = 0. To 
see this, recall that expansion (2) anticipates finite damping and forcing in two ways: it 
assumes the TW modes are in subharmonic resonance with the two driving frequencies, i.e., 
they oscillate at m/2 or n/2, and it takes the wavenumbers k\ and k% from the 7-dependent 
dual minima of the SW neutral stability curves at the bicritical point of the damped forced 
problem not from the dispersion relation, Q,{ki) = m/2 and Q(k^) = n/2. One consequence 
of this implicit 7 dependence in k\ and /c 3 is a detuning in the damped forced problem 
between the parametrically excited frequencies m/2 and n/2 and the "natural" frequencies 
Q(ki) and ^(^3). In Eqs. (8,9) however, the damping and forcing are both zero, and the 
wavenumbers k\ and k 3 coincide with those given by the inviscid dispersion relation, i.e., 
Q(ki) = m/2 and £l(kz) = n/2. Consequently, there is no detuning at / = 7 = 0. 



2.2 Unfolding of Eqs. (8-9): parameter symmetries 



We now unfold the vector field (8,9) assuming all symmetry-breaking parameters are the 
same order: \f m \ ~ \f n \ ~ 7 ~ e. Throughout this unfolding we will be interested only 
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in those terms which add new qualitative effects or determine the leading order scaling of 
coefficients in the subsequent SW equations. Assuming that the symmetry-breaking terms 
are analytic in / and 7, their precise form can be obtained by appealing to a more general 
class of symmetries in which parameters may be transformed along with dynamical variables 
(see, e.g., [22]). In this way the broken symmetries % and a can be recast in terms of 
intact parameter symmetries. Put another way, it is still possible to translate or reflect the 
problem in time so long as the damping and forcing terms are also appropriately transformed. 
Specifically, we have the parameter symmetry 

T s par : Zf) -> (Zf, Z})e ±irns/2 , Zf -> Zfe ±ins/2 , 

(fm,fn)^(fme lmS J n e inS ), (H) 

replacing T s , while a becomes 

o*" : (t, 7 )->-(t, 7 ), Zf^Zj, (f m , f n ) -> (f m , f n ). (12) 

Although parameter symmetries are sometimes unphysical (it makes little sense to take 
7 < 0) they are perfectly good mathematical symmetries of the governing equations and 
may be used to advantage - in this case they yield the proper form of symmetry-breaking 
terms. 



2.2.1 Finite forcing: f ^ 

The presence of forcing allows several new terms in Eqs. (8,9). For example, at linear order 
in amplitude Eq. (8) contains the usual parametric forcing term / m Zf , consistent with the 
symmetry T s par . Another important addition comes in the form of quadratic resonance terms. 
The quadratic amplitude dependence follows from the spatial resonance condition (3) but 
the dependence on f m and f n can only be determined with the help of Z? ar , i.e., it reflects a 
temporal resonance requirement. To determine this dependence we consider adding terms, 
consistent with the spatial resonance requirement, of the form f^f^Z^Zf and f^f^ZfZ^ 
to Eqs. (8) and (9), respectively. Here r ,u, v, and w are integers, allowed to be negative 
to indicate complex conjugation (i.e., f~ 2 = and so on). A combination such as f m f m is 
invariant under T s par and appears only at higher order. The temporal resonance condition, 
which follows from the equivariance requirement under 7^ par , may be expressed as 

m/2 = rm + un — e\m/2 + £2^/2, (13) 
n/2 = v m + wn + e s m/2 + e 4 m/2, (14) 

where £j = ±1, % — 1, ...,4. These distinguish between the Zj modes (the corresponding 
£j are positive) and the Zj modes (negative £j). We seek values for the parameters r, u, v, 
w, £1, £2, £3, and £4 that minimize both |r| + \u\ and \v\ + \w\ and therefore give resonance 
terms of lowest order in /. To accomplish this we rewrite Eqs. (13) and (14): 

m(l + £i -2r) = n(e 2 + 2u), (15) 
m(l + £ 3 £ 4 + 2£ 3 -u) = n(e 3 - 2e 3 w), (16) 
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showing that, as claimed earlier, n must be even else there is no solution at all, and that 
Eq. (16) is equivalent to Eq. (15) under v — * — £2?", w — > —e 2 u, £3 — > £2, and £4 — > £i£ 2 - 
We may therefore focus on Eq. (15), which requires that the quantity 1 + £1 — 2r (hereafter 
A) is a nonzero multiple of n. In terms of A the minimum of |r| + \u\ may be expressed 
(using Eq. (15)) as the minimum of G(A,£i,£ 2 ) = |A — 1 — £i| + \mA/n — e 2 \. For positive 
A this minimum occurs when E\ = £2 = 1 and yields G(A, 1, 1) = (1 + m/n)A — 3 (the 
absolute values can be removed because A > n > 2). Alternatively, if A < the minimum 
of G(A, ei, £2) is realized with e 1 = e 2 = —1 and leads to G(A, —1, —1) = (1 + m/n)\A\ — 1. 
It follows that the global minima occurs for the smallest positive A, i.e., A = n, a choice 
which implies r = — (n — 2)/2 and u = (m — l)/2. The order of the resonance terms (in /) 
is thus \r\ + \u\ = (n + m — 3)/2. 

Altogether, we may rewrite Eqs. (8,9) in the presence of forcing as 

„ _ _ n — 2 m— 1 _ 

Z+ = vZt - (Xf m + zuf^f^Zr + af^f—Z+Zt + 5Z+Z 2 Z^ 
+ (a|Z+| 2 + b|Zf | 2 + c|Z+| 2 + d|Z 2 ~| 2 + e\Z+\ 2 + f|Z 3 -| 2 ) Z+ 

+ (qil^+l 2 + q 2 |^rl 2 + qsl^l 2 + q 4 |^l 2 + q 5 |^ 3 + l 2 + qel^l 2 ) /m^r 

+ {^lfmZt Z 2 + fmi^Zi Z l + ^Z^^ ) + T4f n Z^Z 3 + T 5 f n Z^Z 3 Z^ , (17) 

_ ^ n — 2 _ m — 1 

Zt = QZ+ - (flf n + (fZfn^Zs + Pfm 2 fn 2 Z+Z+ 

+ (g(l^i + l 2 + l^ 2 + l 2 ) + Ml^rl 2 + 1^ 2 1 2 ) + i|^ + l 2 + pI^I 2 ) ^ 
+ ( Sl (|z+i 2 + |z+i 2 ) + S2 (|zri 2 + |z 2 -| 2 ) + S3 |z 3 + i 2 + s 4 |z 3 -| 2 ) / n z 3 - 

+ (ui/ m (Z+Zr + Z+Z 2 -) + u 2 / m (Z+Zr + Z+Z 2 ) + u 3 f n Z+Zz) Z+. (18) 

The remaining equations follow from the symmetries k and 1Z. At second and third order 
in the amplitudes Z^ we have kept all terms which contribute at leading order (in e) to the 
corresponding SW coefficients (Eqs. (32) of Section 2.3). At the linear level, however, we 
have included two additional terms (with coefficients w and g) that do not contribute at 
leading order but are nonetheless responsible for an important qualitative effect, namely, a 
dependence in the linear problem on the relative phase of the two forcing components - this 
dependence will be discussed in more detail below. We note that terms of the form f'Z 3 , for 
example qil^i" \ 2 f m Z^ in Eq. (17), were also found by Zhang [23]. 

Because of time reflection <r par most coefficients in Eqs. (17,18), as in Eqs. (8,9), are imaginary 
at leading order. They do, however, include small /-dependent corrections, invariant under 
>j- par , which can generate real parts. For example, we have 

a = + & m \f m \ 2 + & n \f n \ 2 ) + a/Im(/™/™) + 

where a^, a m , a n , a/ G M.. Analogous expansions apply to all remaining coefficients except v 
and g which contain no 0(1) imaginary part (no detuning at 7 = / = 0). Most /-dependent 
contributions to the coefficients of Eqs. (17,18) are unimportant. In fact, with the exception 
of those in v, g, A and \x (see below) we neglect all of them, keeping only the 0(1) imaginary 
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part and dropping the real parts proportional to Im(f m f™) in anticipation of lower order 
contributions due to finite 7. 



2.2.2 Finite damping: 7 7^ 

The principal effect of the damping 7 is to give real parts to all coefficients in Eqs. (17,18). 
Due to <r par these contributions must be odd under 7 — > —7 and take the form 7, 7 s , 7|/ m | 2 , 
etc. Additional corrections to the imaginary parts, which must be even under 7 — > —7, 
are functions of 7 2 and appear at order e 2 and higher. Such imaginary corrections only 
matter for the detunings Im(fi) and Im(f 3 ), where they contribute at leading order to 
the SW coefficients, and in the (linear) parametric forcing coefficients A and fi where they 
contribute to a (phase-independent) coupling in the SW linear problem between the two 
forcing components f m and f n - some of the 0(e 3 ) terms in the real parts of v and g 
contribute to this coupling as well. Put simply, T W of type Z^ 2 [Zf) are primarily excited 
by f m (f n ) but are affected by /„ (f m ) as well; we are interested in that effect. For the four 
linear coefficients mentioned above we therefore write 

V = 7(-tV + £ 7 7 2 + Vm\fm\ 2 + Vn\fn\ 2 ) + i{v^ 2 + V m \f m \ 2 + V n \f n \ 2 ), V r > 0, (19) 
Q = l(~Qr + £? 7 7 2 + Qm\fm\ 2 + Qn\fn\ 2 ) + KQ^ 2 + Qm\fm\ 2 + Qn\fn\ 2 ), Qr > 0, (20) 

A = A r7 + i{\ + A 7 7 2 + \ m \f m \ 2 + A„|/ n | 2 ), (21) 

jl = /2/y + + /i 7 7 2 + fXm\f m \ 2 + i^n\fn\ 2 ), (22) 

while for the remaining coefficients we use the expansions 

w = ujrj + iwt, a = a r j + ioti, 5 = 5 r ^ + i5 iy etc., (23) 
stopping at lowest order in both real and imaginary parts. 

The form of expansions (19-23) is consistent (where it overlaps) with the analytical formulas 
derived by Topaz and Silber [10] in the case of one-dimensional waves. Those results were 
also for two-frequency forcing in the case of weak damping, and were based on the Zhang- 
Vihals equations. One difference is that Topaz and Silber find no 0(e) real part to A, that 
is, A r = (one would likewise expect fi r — 0). Setting A r = jx r = is of no consequence for 
us, however, because those terms do not contribute at leading order to the SW coefficients 
(see Eqs. (32)). 

It is of interest now to consider the linear problem and, in particular, the effect (mentioned 
above) of the phases (0 m and <p n in Eq. (1)) on the bicritical point defined by 

\v\ = \F m \, F m = A/ m + ^/r7n m , (24) 
|2| = |F B |, F n = jif n + Cf^ff 1 - (25) 

At leading order this bicritical point is independent of <p m and n , i.e., the critical amplitudes 
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|/^J and \fn\ do not couple and we have simply 

\a = ^7, \r n \ = 1,7- (26) 

Expressions (26) are utilized throughout this paper (for example, in evaluating coefficients (19- 
22) at the bicritical point to obtain expressions (A.1-A.7) of the Appendix) and are sufficient 
to obtain the leading order terms in the SW normal form coefficients (see Eqs. (32)). How- 
ever, the lack of phase dependence in Eqs. (26) is contrary to the experimental observations 
of Miiller [12] and the calculations of Zhang and Vihals [15] for the case m:n = 1:2. These 
authors found that there was a strong dependence of r = \fm\/(\fm\ + |/^|) on n . To see 
this dependence Eqs. (24,25) must be solved to high enough order to capture the influence 
of the ^fm^fn an d Cfmfn' 1 terms. When this is done for mm = 1:2 the expression for r 
is 

-7COS0/, (27) 



\2 



V r \fri\ + Qr\K\ |Ai|(lv|/Zi| + Qr\K\y 

which reveals an 0(e) harmonic dependence on the relative phase 

0/ = m(f) n - n<p m . (28) 

This particular linear combination of m and <f> n , the complex phase of the 7^ par -invariant 
term /™/™, is the only physically relevant one. In general, the ^/-dependent correction to 
r is (9( e m +™- 2 ) a nd is much less important (with small damping) when m + n is large. 
For the next available resonant triad interaction, m:n = 3:2, it is already of third order. 
The diminishing influence of 0/ with m + n is consistent with the experiments of Edwards 
and Fauve [4] where no significant phase dependence was found for m:n = 5:4. In the case 
m:n = 1:2 we can compare Eq. (27) with [15], where a different convention for the two- 
frequency forcing function is used, by setting <p rn = tt/2 and n = n/2 + <fi (0 is the phase 
variable used in [15]) to find <pf = — 7r/2. Thus Eq. (27) predicts an extremum for r(0) at 
= 7r/2, in agreement with the results of [12] and [15] (see Fig. 1 of [15]) - to be more precise, 
a maximum was observed at = ir/2. That a maximum should be observed rather than a 
minimum cannot be determined from Eq. (27) since this depends on the actual numerical 
values of the coefficients v r , fii, etc.. 



2.3 Reduction to standing wave equations 



The 0(e) linear coupling between and ZJ due to the parametric forcing destroys the pure 
TW solutions of Eqs. (8,9) and ensures that the first instability with increasing / produces 
SW satisfying \Z~j~\ = \Z~\. A standard reduction procedure (see, for example, [24]) may 
therefore be used to reduce Eqs. (17,18), near the bicritical point, to a system of three 
amplitude equations describing SW modes. Small deviations from the bicritical point are 
captured with the reduced forcing parameters A m , A n e R according to f m = /^(l + A m ), 
fn = /^(l + A n ). To facilitate the reduction we first rescale the amplitudes so that (after 
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appropriate redefinition) the complex phases of F m and F n are equal, respectively, to those 
of v and g. Specifically, this rescaling takes {Zf, Z 2 ) — > (Zf , Zf) exp{±i$i/2} and Z 3 — > 
Zf exp{±i$ 3 /2}, where -F^t) = \F^ n v c \exp{i r &\\ and i^£ c = \F^g c \exp{i'd 3 }. In addition 
to F rn and F n , phase shifts are induced in all of the nonlinear coefficients in Eqs. (17,18) 
except a-h, 1, and p. As a result of this rescaling the critical eigenmodes satisfy the simple 
relation Zf = Z~ = Aj. 

After the reduction we find that the Aj satisfy 

A 1 = \Ai + aA 2 A 3 + A l (a\A^ 2 + b\A 2 \ 2 + c\A 3 \ 2 ) , (29) 
A 2 = \A 2 + aA,A 3 + A 2 (a\A 2 \ 2 + b\A^ 2 + c\A 3 \ 2 ) , (30) 
A 3 = fxA 3 + j3AiA 2 + A 3 (rf|A| 2 + d\A 2 \ 2 + e\A 3 \ 2 ) , (31) 

where 



A 



v r jA m + J l7 3 A n + e x wi (mA n + (n - l)A m ) \ f^'^fX cos < 



/i = g r >yA n + J 2 7 3 A m + e^Q (nA r 



[m 



a = -e x ai\f^\ n 2 |/, 

~ -n —2 

a = J37, 



m— 1 



m— 1 
2 



COS $, 
COS $, 



6 = J 4 7 - 
c = J57 + 

d = J&l - 
e = J 7 7. 



2^7 

2^r7 

2iv 7 



1 f° \ n ~ 
\J m\ 


- 2 \fT~ 


1 sin 2 $, 


\ rc \n- 
\J m\ 


- 2 \fT~ 


1 sin 2 $, 


1 rc \n- 
\J m \ 


- 2 \fT- 


_1 sin 2 $. 



(32) 



Here Ji, J 7 are 0(1) functions of the TW coefficients (19-23) and are given in the Ap- 
pendix, e x = sign(Aj), = sign(/2j), and 

* = £ «J - f ■ (33) 

The coefficients (32) are all real, expressed in terms of the unsealed parameters (19-23) 
and include the lowest order in e. Additional interesting contributions, not necessarily at 
lowest order, have been kept as well. In A and fi, for example, we have included the leading 
order nonresonant (i.e., lacking explicit dependence on m, n, and 0/) coupling between 
the two forcings (the Ji7 3 A n term in A and the J 2 ^A m term in /i) and the lowest order 
"resonant" contribution (the (9(e m+n_1 ) terms varying as cos0/). In the nonlinear cross- 
coupling coefficients b, c, and d there is an 0( e m + ri - 4 ) resonant contribution coming from 
the damped SW modes that have been slaved away; this contribution appears at leading 
order when m + n < 5 but should be negligible for m + n > 5. 
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We now summarize the most important properties of the SW coefficients (32). All of this 
structure relies on the parameter symmetries 7^ par and <r par , both possessed by the full 
hydrodynamical problem [25], the restriction to weak damping (7 ~ \f^\ ~ \f£\ <C 1) and 
the assumed analytic parameter dependence on 7 and /. 

(1) For m + n > 5 all coefficients are 0(e) or higher. A similar result was 
found in the analytical expressions of Topaz and Silber [10] and of Zhang 
and Vinals [15,23]. The latter authors point out [14,23] that the appearance 
of 7 in the nonlinear coefficients despite the fact that the Zhang- Vinals 
equations contain only linear dissipation, should not be surprising. There 
is no reason why parameters in the linear terms of the governing equations 
cannot be involved in the nonlinear coefficients of the amplitude equations. 
In our derivation this dependence is a natural consequence of the underlying 
temporal symmetries. From a physical perspective it is also not surprising 
that the timescale for dissipative dynamics diverges as I/7. 

(2) The unfolding parameters A and fi are linear combinations of A m and A n , 
and are therefore rotated slightly (and stretched) with respect to the latter, 
"more natural" parameters. This rotation is 0(e 2 ) unless m:n = 1:2 in 
which case it is 0(e) and proportional to cos 0/. Similar rotation is discussed 
in more detail in [10], although those results are obtained away from the 
bicritical point. 

(3) When m:n = 1:2 the cubic cross-coupling coefficients b, c, and d depend (si- 
nusoidally) on the relative phase 0/ at leading order. This behavior is similar 
to the leading order phase dependence found by Zhang and Vinals [15] for 
cubic cross-coupling coefficients with 1:2 forcing, a result obtained perturba- 
tively and not from symmetry considerations. Furthermore, the coefficients 
b, c, and d diverge as e — > 0, indicating a failure of the reduction procedure 
in this limit. This failure is not surprising because at 7 = there are three 
additional critical modes which need to be included. The three neglected 
SW modes satisfy Z^ ~ — Zj after the rescaling prior to Eqs. (29-31) and 
experience 0(e) damping. The exceptionally strong coupling to these other 
modes when m:n = 1:2 (the resonance terms of Eqs. (17,18), through which 
the coupling occurs, are 0(e) or smaller unless m:n = 1:2) helps explain 
the failure of the reduction procedure leading to Eqs. (29-31). In general, 
the description based on three SW modes is not expected to be valid if 7 is 
"too" small. 

(4) All measurable effects from coefficients (32) are 27r-periodic in 0j. This is 
true even for a and f3 (which appear to be 47r-periodic in 0j; see Eq. (33)) be- 
cause, owing to the invariance of Eqs. (29-31) under (A,-, a, (5) — > — (A,-, a, (3), 
only the relative sign of a and j3 matters; the dynamics with $ and $ ± n 
are equivalent. 

(5) The resonant coefficients a and (3 are of order e ( m + n - 3 )/ 2 anc [ thus depend 
strongly on the choice of forcing frequencies, their influence diminishing 
exponentially with m + n; this fact helps explain why this type of resonant 
triad is not seen for large values of m + n [3,8]. The e ( m + n - 3 )/ 2 dependence is 
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not universal however. Near particular choices of <p m and <p n where cos <£> = 
the terms given in Eqs. (32) vanish and, assuming the next nonvanishing 
term (not calculated) is a factor of 7 smaller, there is a cross-over to the more 
extreme e ( n + m - 1 )/ 2 dependence. Conversely, one could say that resonance 
effects are strongest when $ = 0, ir. 



2.3.1 Consequences of Hamiltonian structure 

In this section we examine the consequences of full Hamiltonian structure (not just the time 
reversal symmetry <r par ) in the undamped (i.e., inviscid) problem; see [14,17,26,27,28,29,30]. 
Specifically, we suppose that Eqs. (17,18) possess a Hamiltonian H with the complex am- 
plitudes Z^ obeying the equations of motion (see, e.g., [31,32]) 

Zf^Ti-^E, J = 1,2,3. (34) 

Requiring that Ti be a real scalar function, invariant under the symmetries 7^, 1Z, k, Tf w , 
and a par leads to 

n = hdzfi 2 + \z±\ 2 ) + h 2 \zf\ 2 + h 3 (f m (z+Zi + z+z 2 ) + c.c.) 

+ h 4 (f n z+Zz + c.c.) + ^(/rvr^r + Z 2Z 2 ) + c.c.) + h % {rjr x ztz 3 + c.c. 

_ n — 2 m — 1 _ _ n — 2 _ m — 1 _ _ 

+ h 7 (frn 2 fn 2 Z+ 2+ Z+ + f m 2 fn 2 Zf Z 2 Z3 + C.C.) + h 8 (\Z+ \ 2 \Z[ | 2 + \Z+\ 2 \Z 2 \ 2 ) 

+ h 9 (\Zt\ 4 + \Zt\ A ) + h w \Zf\ 2 \Zf\ 2 + hn\Z?\ 2 \Z*\ 2 + h l2 \Zf\ 2 {\Zf\ 2 + \Z±\ 2 ) 
+ h^ZfZ^Z^Z^ + /ii 4 |^| 4 + /nsl^l^lZfl 2 + \Z^\ 2 ) + /i 16 |Z 3 + | 2 |^| 2 

+ h 17 (f m ZtZ 2 -{\Zt\ 2 + |Zr| 2 ) + / m Z+Zr(|Z 2 + | 2 + \Z 2 \ 2 ) + c.c.) 

+ h 1B (f m Z+Zr(\Z+\ 2 + |Zf| 2 ) + / m Z 2 + Z 2 -(|Z+| 2 + |Z 2 -| 2 ) + c.c.) 

+ h 19 (f m (Z+Z- + Z 2 +Z 2 -)(|Z+| 2 + |Z-| 2 ) + c.c.) + fc 2 o(/ n Z 3 - Z+(|Z 3 + | 2 + \Z~\ 2 ) + c.c 

+ /i 21 (/ n Z 3 + Z 3 -(|Z+| 2 + |Zr| 2 + |Z+| 2 + |Z 2 -| 2 ) + c.c). (35) 

Only those terms needed for comparison with Eqs. (17,18) have been included in Eq. (35). 
The coefficients hi, h 2 \ are real, and those terms involving a superscript ± should be un- 
derstood as a sum over both signs, e.g., \Zf\f = \Zf | 2 + |Zf | 2 . The equations of motion (34) 
are equivalent to Eqs. (17,18) only if 

a = (3, e = g, h = f, qi = 2q 2 = 2r 2 , q 3 = q 4 = n = r 3 , 

Q5 = Qe = ui = u 2 , r 4 = r 5 = si = s 2 , s 3 = 2s 4 = 2u 3 , (36) 

where all coefficients are evaluated at 7 = 0. The most important of relations (36) is a = $ 
(i.e., oti = /3i) because that implies that the resonant coefficients a and (3 of Eqs. (32) have 
opposite signs for small 7. Although this situation can change with increasing damping, the 
fact that there is good reason to expect a(5 < has profound dynamical implications for 
Eqs. (29-31) - it is a prerequisite for a variety of interesting drifting solutions and heteroclinic 
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cycles (see Section 4) . Chossat [33] has shown that the special form of advective nonlinearities 
leads to the same conclusion (resonant coefficients of opposite sign) for a general class of 
self-adjoint hydrodynamical problems (see also [34]). A strict interpretation of the equation 
aij = Pi (which would also conclude that a and (3 have equal magnitudes; see Eqs. (32)) 
is inappropriate here because, for one thing, such a statement depends on normalization 
conventions. Provided a and (3 do not vanish, the amplitudes Aj may always be rescaled 
to set \a\ = \f3\ = 1. More importantly, to say that the TW equations (17,18) display 
"Hamiltonian structure" need not imply that the Zf and Z^ are themselves canonically 
conjugate Hamiltonian variables in the sense of Eqs. (34). Transformations (simple rescalings, 
for example) that preserve dynamics should be allowed as well. With this relaxed criteria 
we predict for the cubic coefficients simply that the oscillations of b, c, and d with $ are 
the same order in magnitude and that the oscillating part of b and of d (measured with 
respect to the nonresonant 0(e) part) is strictly negative while that of c is positive. These 
oscillations appear at leading order only if m:n = 1:2 or 3:2. In the former case they are 
0(1/ e) and can overwhelm the $- independent 0(e) part. 



2.3.2 Discussion of standing wave equations (29-31) 

We have related many properties of the coefficients of Eqs. (29-31) to the influence of weakly 
broken temporal symmetries and Hamiltonian structure. The spatial symmetries of the prob- 
lem, however, remain fully intact and must be considered as well. They act on the SW 
(cf. Eqs. (3-5)) according to 

T : Aj -> A je ^, = (0 l5 2 ), <j> u 2 G [0, 2tt), 3 = 0i + 02, (37) 
k : A 1 <-> A 2 , (38) 
H:Aj-> Aj, (39) 

In particular, the action of 1Z explains why the coefficients in Eqs. (29-31) are real. Note that 
this nontrivial action of 1Z on the SW amplitudes, a necessary consequence of the spatial 
symmetries of the problem, allows for the possibility of a symmetry-breaking bifurcation to 
traveling (hereafter, drifting) waves. Strictly speaking, the "standing waves" of the previous 
section are standing only on the fast timescale set by u (i.e., the eigenfunctions satisfy 
\Zj~\ = \Z~\). On the slow timescale they may or may not drift, depending on whether or 
not they belong to the 7£-invariant subspace (or a translation of it). 

Recall that unless n is even there are no quadratic resonant terms, i.e., a — f3 — 0. This 
result, obtained in [11], follows directly from Eqs. (29-31) by noting that the forced problem 
is invariant under a discrete time translation through one period of the forcing - this is the 
extent to which time translation is usually considered. Accordingly, Eqs. (29-31) must be 
equivariant under the operation (A 1: A 2 ) — > (— l) m (A 1 , A 2 ), A 3 — > (— l) n A 3 . If n is even (and 
therefore m odd) this discrete temporal symmetry is equivalent to the spatial translation 
T{k^) and imposes no additional restrictions [11]. However, if n is odd the induced symmetry 
forces a — (3 — 0. Because we are interested in resonance effects we always take n to be 
even, i.e., we assume A 3 is harmonic with respect to the overall forcing period 2n. 
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Even when A3 is harmonic and the resonant coefficients are nonzero there is no guarantee 
that resonance effects can be easily observed. As seen in Eqs. (32), this depends strongly on 
m + n and 7. If a and (3 are extremely small, their influence will be appreciable only while 
the system is extremely close to onset. Although in one sense this limit is ideally-suited 
for the weakly nonlinear approach taken here, it may require far too much sensitivity for 
realistic experiments, i.e., the smallest available turn of the knob may immediately put the 
system into a regime dominated by cubic nonlinearities. To estimate the likelihood of this 
occurring we characterize the range over which resonant terms contribute relative to the 
range over which the weakly nonlinear analysis is likely to be valid. We may assume, for 
example, that a weakly nonlinear approach is "justified" when \Aj\ < e , where e = 0.1 
(or any other reasonable value). The dynamics, however, will begin to switch from resonant 
(quadratic nonlinearities) to nonresonant (cubic nonlinearities) when these two effects are 
of the same order: (a, (3)\Aj\ 2 ~ (a, b, c, d, e)|A,| 3 . For concreteness, we suppose that this 
transition occurs when \Aj\ ~ e rcs = max(|a|, |/3|)/max(|a|, \b\, |c|, \d\, |e|) and denote the 
ratio of e res to e by r] res : 



max ( a 


.1/31) 


eo max ( a 


,\b\, 


c 


,141 


e|) 



(40) 



Thus, resonance effects are expected to be important throughout the weakly nonlinear regime 
if r) res > 1, while an earlier crossover to nonresonant behavior is anticipated when ?7 res < 1. 
Note that unless m:n = 1:2 (this case is more complicated because b, c, and d diverge as 
e — > 0) r] ICS scales as e O+ n - 5 )/ 2 . Broadly speaking, for small 7 (i.e., e) one expects that 
resonance effects will be easily observed when m + n < 5 but will become increasingly more 
delicate as m + n increases. 



3 Numerical results using Zhang- Vinals equations 

In this section we test the predictions of Eqs. (32) by numerically calculating the SW co- 
efficients from the Zhang- Vinals model [14,15], a quasipotential formulation of the Faraday 
problem capturing the small-amplitude dynamics of deep fluid layers in the limit of weak 
damping. This calculation also gives physically relevant coefficients which we use in Sec- 
tion 4 to illustrate the SW dynamics. We emphasize that the procedure described here is a 
direct reduction of the Zhang- Vinals equations at the bicritical point to three SW amplitude 
equations, and does not rely on the symmetries of Section 2.2. 

In the frame of the vibrated fluid we can write the effective gravitational acceleration in the 
form g(t) = g — g m cos(mcjt + m ) — g n cos(nut + <p n ) , where g > is the usual gravitational 
acceleration. Note that we are picking a convention by using cosines rather than sines, and by 
putting an overall minus sign in front of the applied acceleration. Switching to the opposite 
sign convention is equivalent to taking <p m — > <p m + it, <p n — > <p n + it. Because of Eq. (28), we 
then have 0/ — * 0/ + (m — 7i)tt and $ — > $ — (m — n)ir jl. Since most of the normal form 
coefficients in Eqs. (32) depend on <pf this convention does matter. In a, f3, b, c, and d it 
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amounts to a shift by it/2 in the oscillations versus This follows from the fact that m — n 
is odd, together with the invariance of the problem under $ — > $ ± 7r (i.e., 0/ — > 0/ =F 27r; 
see Section 2.3); changing the above sign convection is thus equivalent (for a, (3, b, c, and 
d) to taking $ — > $ + n/2. The difference is observable - in fact, many other patterns, e.g., 
hexagons, depend on a phase analogous to $ (see [9]) - but, to our knowledge, even those 
experiments [3,4,6,8,12] which investigate dependence on the temporal phase (equivalent to 
(f)f) fail to report the convention used (one may infer the convention used in [12] from the 
agreement found in [15]). Without this information, such results are ambiguous. 

The governing equations describing surface height h(x,r) and surface velocity potential 
ip(x,r) are expressed in nondimensional form as 

d T h = 7 V 2 /i + V(p-V- (hWip) + -V 2 {h 2 V v ) ~ T>(hV(p) + v(hV(hV(p) + -/i 2 VV), 
d rV = 7 VV + ToV 2 /* - G(r)h + l -(V^Y - i(Vy?) 2 

-(%)(wV + «) -^r v- ((v/*)(v/i) 2 ), (4i) 

where G(t) = G - f m cos(mr + m ) - f„ cos(nr + n ), and (f m , f n ) = (^ m , g n )k /u 2 ; 7, G 
and T are as defined in Section 2. The nonlocal operator V (see [14]) multiplies each Fourier 
component by its wavenumber: Vexp{ik ■ x} = kcxp{ik ■ x}. 

To extract coefficients from Eqs. (41) we use a three-timing perturbation technique (see [16] 
or [11] for more details), introducing a bookkeeping parameter e and writing 

h = ehi + e 2 h 2 + ^ = e<pi + e 2 Lp 2 + (42) 

f m = f^ + e 2 f m , f n = f c n + e 2 ~f n , d T ->d T + ed Tl +e 2 d T2 + .... (43) 

After hi and ip 1 are expressed in terms of Fourier modes, e.g., /i 1 (x,r) = J2kPk( T ) e w{^' x } 
and </?i(x, r) = J2k Q , fc( r )exp{ifc ■ x}, the 0(e) problem reduces to a series of uncoupled 
damped Mathieu equations parameterized by the wavenumber k: 

d TT p k + 2'jk 2 d T p k + (7 2 /c 4 + G k + T k 3 )p k = fc(f m cos(mr + <p m ) + f n cos(nr + <p n ))p k 

(44) 

Note, from the first of Eqs. (41), that q k = k~ x (d T + ^yk 2 )pk- At the bicritical point (f m , f n ) = 
(f^,f°), there are exactly two periodic solutions to Eq. (44), p^ij) and pk 3 (r), associated 
with the dual minima of the neutral stability curves (see Fig. 2). 

It is worthwhile at this point to extract the leading order linear TW coefficients in Eqs. (17,18) 
directly from Eq. (44). First, by setting 7 = and treating f m and f n as small perturba- 
tions, Eq. (44) can be solved with a multi-timing scheme that gives, as the first solvability 
condition, 

Zf = -ik 1 /{2m)l c m Zf, Zt = -ik 3 /(2n)i c n Z*. (45) 
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Upon comparing Eqs. (45) with Eqs. (17,18) and (23) we identify 



A,; = 



2m'' 



IH = 



2n 



(46) 



A complementary treatment with f m = f n = and 7 in the role of small perturbation 
establishes, via Eqs. (17,18) and (19,20), that 



k 2 



Qr — ^3) 



(47) 



Eqs. (46) and (47) agree with the formulas given in [10]. In Eqs. (47) we see that the linear 
damping is proportional to k 2 , as expected. The significance of Eqs. (46), on the other hand, 
is in relation to the forcing convention discussed above. We have found that Aj > and 
jli > 0, i.e., that E\ — — 1 in Eqs. (32,33). Had the alternative convention been used 
where g(t) = go + g m cos(mut + <f> m ) + g n cos(nut + <f> n ), we would have e\ = = — 1. 

Since we are assuming n is even and m is odd, p kl is subharmonic with respect to the forcing 
period 2n: p kl (r + 2n) = —p kl (r), while p k3 is harmonic: p kz (j + 2n) = p k3 (r). To focus on 
resonant triads, we expand hi and ipi in terms of the three spatially resonant critical modes 
(see Fig. 1): 



hi = p k ,{r)[Ai{T u T 2 )e ik ^ + A 2 (T 1 , T 2 )e ik2 ' x 
ipi = q kl {T)\Ai{Ti,T 2 y k ^ x + A 2 (Ti,T 2 )e ik2X 



+ C.C. 



+ p k3 {r)A 3 (Ti,T 2 )e 
+ q k3 (r)A 3 (T u T 2 )e^- x + c.c., 



(48) 
(49) 



At 0(s 2 ) and 0(e 3 ) one obtains differential equations from the solvability conditions. These 
can be reconstituted by dropping e and setting d t = <9ti + <9t 2 - The subsequent equations 
are identical in form to Eqs. (29-31), as they must be from symmetries (37-38). The ex- 
plicit expressions for these coefficients are in general quite lengthy and must be computed 
numerically since they involve Fourier representations of p k (r) and q k (r). For this reason we 
present only numerical results, concentrating mainly on the case m:n = 3:2 for which there 
are many recent experimental results [3,5]. 

When m:n = 3:2 the coefficients a, /3, b, c, and d depend at leading order on $. We find, 
however, that for the Zhang- Vihals equations the ^-independent parts of b, c, and d are 
considerably larger than the ^-dependent parts despite the fact that they formally appear 
at the same order in 7. The most important effect of $ is therefore in the cos $ behavior of 
a and (3 (see Eqs. (32)) which we confirm in Fig. 3. Observe that «($) and /3(<&) have the 
opposite phase, i.e., af3 < 0, as predicted by the arguments of Section 2.3.1. Maximal values 
of I a I and \f3\ are realized (for small 7) when $ = 0. We use such an "optimal" choice of 
phases, e.g., m = — tc/4 and 4> n = 0, to test in Fig. 4 the scalings predicted by Eqs. (32). 
Note that, as expected for m:n = 3:2, we find r] ICS ~ 0.02-0.3, and resonance effects are 
likely to be important over much of the weakly nonlinear regime. 

For values of $ near ±7r/2, the 7 dependence of a and (5 switches from the exponent 
(m + n — 3)/2, characteristic of the terms given in Eqs. (32), to the exponent (n + m — l)/2, 
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-7T -IT/2 7T/2 7T -7T -7r/2 7r/2 7T 



Fig. 3. Dependence of a and /3 on <3? in the case m:re = 3:2 with physical parameters 
(in cgs units): p = 1, v = 0.1, T = 21.5, and u = 10, corresponding to 7 ~ 2.08x10 . 
To change <3? we set 4> n = and varied </> m . 




i(T 5 icr 3 io~ 1 icr 5 -icr 3 icr 1 icr 5 icr 3 icr 1 



Fig. 4. Behavior of 7(0;), r] res (cu) (see Eq. (40)), and the magnitude of the coeffi- 
cients | a (7) |, | ^(7) | of Eqs. (29-31) over the interval u £ [3,300] rad/s with 
4> m = — vr/4, (j) n = 0, p = I, v = 0.1, and T = 21.5 (cgs units). The straight lines 
show the predicted (linear) 7 dependence. In the calculation of rj res we use eo = 0.1. 

reflecting the next nonvanishing term (not calculated). Fig. 5 illustrates this second type of 
behavior by showing T] ieB {co), |a(7)|, and |/3(7)| when $ = ir/2. The remaining coefficients 
in Eqs. (29-31) are essentially unchanged from Fig. 4. Note that r] res — > as u — > (i.e., 
7 — > 0), indicating that resonance effects are becoming increasingly weak - larger values 
of u are required to see them more easily. Furthermore, besides the fact that a and (3 are 
smaller, we find that a/3 > over the calculated range ui G [3, 300]. This is in sharp contrast 
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Fig. 5. Dependence of rj Tes (uj) (using eo = 0.1), |a(7)|, and \{3(i/)\ when <E> = tv/2. 
Parameters are the same as in Fig. 4 except with <p m = tt/4. Straight lines show the 
predicted 7 2 dependence. 

to the typical case, with $ bounded away from ±7r/2, when the leading order terms given 
in Eqs. (32) dominate and af3 < 0. Because the relative phase 0/ (equivalently, $) can alter 
the sizes of a and (5 and the sign of their product, it can have a dramatic effect on the 
dynamics of resonant triads, as shown in the following section. 



4 Dynamics 

In this section we discuss the dynamics of Eqs. (29-31) with an eye to the influence of 
relations (32) and (36). As shown by Guckenheimer and Maholov [35] these dynamics can 
be very rich, in part because a number of solutions drift or "travel" on the slow time scale, 
as discussed in Section 2.3.2. In the following we refer to these as drifting waves (DW) to 
emphasize that this motion is much slower than the fast (standing) oscillations associated 
with the parametric forcing. 

While the temporal symmetries are now implicit, expressed in relations (32) and (36) de- 
scribing the scaling, parameter-dependence, and relative sign of the SW normal form coeffi- 
cients, the unbroken spatial symmetries (37-38) force the existence of dynamically invariant 
subspaces containing distinct classes of solutions; the most important of these fixed point 
subspaces are listed in Table 1. Not all of these subspaces are isolated and independent. For 



Subspace 


Fixed by 


Defining property 




K 


Im(A,-) = 


u K 


K 


A 1 = A 2 


Ui 


T (O,0) 


A 2 = A 3 = 


u 2 


%0) 


A 1 = A 3 = 


Us 




A 1 = A 2 = 



Table 1 

Invariant subspaces of Eqs. (29-31). 

example, Un contains the solutions that are invariant under inversion through the origin 
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(i.e., under 1Z) but there is a torus of equivalent subspaces generated by 7^; each of these 
related subspaces is invariant under inversion through an appropriate (translated) origin. 
Similarly, there is a circle of subspaces equivalent to IA K generated by solutions in 

these subspaces satisfy A 2 = e l2<t> Ai and are reflection-symmetric about a line parallel to k 3 . 
The subspaces U\ and U 2 are related by the symmetry k and are equivalent. In the following 
we use Uft, U K , and U\ to denote the entire family of equivalent subspaces, the representative 
members of which are listed in Table 1. In practice, the most important subspaces are Un 
and U K because they relate to easily recognizable characteristics: solutions in Un are strictly 
standing (i.e., they do not drift; see below) while solutions in U K possess a line of symmetry 
and may drift only in a direction parallel (or antiparallel) to k 3 . 

For most purposes it helps to factor out the continuous translation symmetry by writing 
Aj = cije l9j and introducing the 7^-invariant phase 9 = 6> 3 — 9± — 9 2 . One then obtains the 
four-dimensional system 



f\ = Ari + ar 2 r 3 cos 9 + ri(ar\ + br\ + cr\), (50) 

r 2 = Ar 2 + ar\r 3 cos9 + r 2 {ar\ + br\ + cr^), (51) 

f 3 = fir 3 + fir\r 2 cos 9 + r 3 (drl + dr\ + er^), (52) 
A (o r \ r 2 . r 2 r 3 rxr 3 N . 

9 = — [p ha ha sin 6^, (53) 

v r 3 ri r 2 / 

where the 9j satisfy 

d r2r3 • a a rir3 • a a n r ^ 2 - a (ka\ 
9i = a smB, 9 2 = a sm9, 9 3 = —p sm9. (54) 

n r 2 r 3 

Eqs. (54) guarantee that solutions in tin, which corresponds to {9 = 0} U {9 = n} in 



Eqs. (50-53), are also standing (i.e., all phase velocities are zero). Furthermore, from Eq. (53) 
it follows that unless a(3 < the quantity multiplying sin 9 is single-signed and all trajectories 
approach U-ji (see [35,36]), i.e., there are no permanently drifting solutions if a(3 > 0. The 
fact that one can expect drifting solutions on the basis of underlying Hamiltonian structure 
was the conclusion of Section 2.3.1 

4-1 Principal solutions 

We describe here the most important solutions of Eqs. (29-31), equivalently, Eqs. (50-52). 
Many of these solutions (or their analogs) are discussed in [35] . They are summarized, along 
with the bifurcations they undergo, in Table 2. 

(1) Standing rolls SRi G U x 

The SRi states satisfy |^4i| 2 = —A/a, A 2 = A 3 = 0. Physically, they represent SW with dom- 
inant frequency m/2 described by a single wavevector k x . Standing rolls SR 2 with wavevector 
k 2 are related to SRi through the symmetry k and are equivalent. Both bifurcate from the 
flat state (O) when A = 0, supercritically if a < 0. The eigenvalues of SRi within the 
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u K 


u K nu n 




,T'--DW« — - 

* HET2 


-"MM K v. 
MS\\„: 
~~~HETi 


MM-''V; 

MSW''/ 



no symmetry 



DW 



Table 2 

Typical solutions (see text) classified according to the reflection symmetries 1Z and 
k. Broken lines indicate potential bifurcations. 

two-dimensional subspace U\ are (a result of the translational symmetry T^) and — 2A. 
Transverse to U\ the eigenvalues £ have multiplicity two and satisfy 



a 2 £ 2 - a (a/ji + A(a - b - d)) £ + A(a - 6)(a// - dA) + a/5A = 0. 



(55) 



Near onset (i.e., |A| <C 1) these eigenvalues are approximated by /i and A(a — 6 + a(3/fi)/a. 
Thus, the standing rolls SRi are stable at onset if a < (the bifurcation is supercritical), 
fj, < (stability with respect to perturbations of type A 3 is inherited from the basic flat 
state), and /i(b — a) > a(3; they are unstable otherwise. 

The SRi states undergo a steady-state bifurcation along the line 



ad d , 
+ -A, 



b — a a 



(56) 



giving rise to a branch of mixed modes (MM) contained in U-ji with A 2 7^ 0, A 3 7^ 0, and 
experience a Hopf bifurcation along the line 



H = 



b + d 



A, 



(57) 



provided af3 < X(a — b) 2 /a < 0. This Hopf bifurcation with 0(2) symmetry (see, e.g., [18])) 
simultaneously produces modulated standing waves (MSW) contained in Un, and DW not in 
Un- Note that a/3 < is a necessary condition for this bifurcation, as it must be if DW are 
produced. Assuming this condition holds, the line of Hopf bifurcations (57) intersects the line 
of steady bifurcations (56) at the Takens-Bogdanov point: (A, jj) = afl(a, b + d — a)/(b — a) 2 . 

(2) Standing Rolls SR 3 G U 3 . 



The SR3 states satisfy Ai = A 2 = 0, | 1 2 = — /i/e and represent harmonic SW (dominant 
frequency n/2) with wavevector k 3 . They bifurcate from O when \i = 0, supercritically if 
e < 0. Within U 3 the eigenvalues are and — 2/x while the remaining four eigenvalues come 
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in pairs given by 




(58) 

In contrast to SRi, only steady-state bifurcations can occur on SR3 and these generate 
branches of mixed modes MM K G U K (i.e., the unstable eigenvectors satisfy |Al| = \ A 2 \). 



(3) K-symmetric mixed modes MM K G Un fl U K . 

These mixed modes bifurcate, along with SRi, from O when A = 0. The 7^-invariant phase 
9 must be or it while the amplitudes satisfy 

(2cd - (a + b)e)rl ± (2ad + (3c)rl + (2d\ + a(3-(a + b)n)r 3 ± (5\ = 0, 

2 2 cr^iarg + A 

r x = r 2 = * — , (59) 

a + 

with cos# dictating the ± sign. Eqs. (59) show that there can be up to three distinct MM K 
states for each of the two possible 9 values. The additional solutions arise in saddle-node 
bifurcations, and other transitions can occur as well. Hopf bifurcations lead to modulated 
standing waves, MSW K , still contained in U K fl U-r, while symmetry-breaking bifurcations 
lead to drifting waves DW K € 1A K (broken TZ symmetry) or to mixed modes MM e U-n (broken 
k symmetry). The MM states bifurcate only from the MM K branch with cos 9 = sign(o;(a— b)) 
and do so when 

- + -A + y (aa ± Ja^-\c(b-a)A. (60) 
a c c z (a — bY v v 7 



a(3 

V = - 



(4) /t-symmetric drifting waves DW K G U K . 

These solutions require af3 < 0. Within Eqs. (29-31) they drift with constant velocity parallel 
to ±fc 3 (i.e., 9\ = 9 2 = O3/2). Within the four-dimensional reduced system (50-53) they 
appear as fixed points with 9 7^ or n satisfying 



T3 Aa(a + b + d)-f3(2c + eY Tl T * (3 ' 3 ' 

r ^ cos0 - a {Aa{a + b + d)-(3{2c + e)) ' (61) 

and exist in the region 

(/i(2a(a + b) - c{3) - X(4da - e{3)) 2 < (fj, + 2\)a 2 (4a[3(a + b + d)- (3 2 (2c + e)). (62) 

The DW K states themselves can undergo Hopf bifurcations to create modulated drifting 
waves MDW K G U R , as well as symmetry-breaking bifurcations to DW not in U K . This latter 



,2 



23 



bifurcation occurs along the line 

4a(b + d-a) -fie 

f i = A Q ^ ( 63 ) 

Aaa - (3c 

(5) Mixed modes MM G U n 

These mixed states (without k symmetry) satisfy 

9 a(3a + (au — d\)(a — b) „ . . . ,., 

rl= , -IT?-, \ cos0 = sign a a-fc , 

(a — o)(ca — ae) 



2,2 A + crf 

r l + r 2 = » r l r 2 = 

a 



CK 



r 3 . (64) 



They form a branch of solutions connecting the SRi (SR2) states with the MM K branch 
having cos# = sign(a(a — b)). There are no saddle-node bifurcations but if a(3 < there is 
an 1Z symmetry-breaking instability that creates drifting waves DW along the line 

e. ad(bc + cd — ac — ae) 
»=c X + cJa-^W • (65) 

(6) Drifting waves DW. 

These steadily drifting solutions (requiring a(3 < 0) of Eqs. (29-31) are not in U K and 
are therefore free to move in any direction (i.e., their drift direction is determined by the 
coefficients, not by symmetry). They are represented by fixed points within Eqs. (50-53) 
satisfying 

^3 _ op + (a - b - d)\ ji , ji _ A + erg 
bd + cd — ac — ae a 

2/j (a -6) 2 (A + erf) ar 3 cos^ 

cos 6 1 = , r 1 r 2 = — . (66) 

apa a — b 

As with MM (and DW K ) there are no saddle-node bifurcations. In fact, for the sets of pa- 
rameters we investigated there were no additional bifurcations on the DW branch. 



(7) /t-symmetric modulated standing waves MSW K G Un D IA K . 

These states appear via Hopf bifurcation from MM K . The form of the bifurcation set is 
nontrivial but reduces to a parabola (see [37,38]) near the origin (//, A) = (0, 0): 

3e 

H ~ -— rA 2 , a(3e\ < (67) 
\a\ 

The MSW K solutions may undergo an 1Z symmetry-breaking bifurcation to modulated drift- 
ing waves MDW K and are typically destroyed in a global bifurcation involving a heteroclinic 
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connection between O and SR 3 (see [36,39]). 

(8) Modulated standing waves MSW G Un . 

These solutions do not have the k symmetry of MSW K . They appear, along with DW, in a 
Hopf bifurcation on SRi. 

(9) /t-symmetric modulated drifting waves MDW K G U K . 

These three-frequency states can emerge from DW K through Hopf bifurcation or from MSW K 
via an 1Z symmetry-breaking bifurcation. 

(10) Heteroclinic cycles HET. 

There are (at least) three kinds of heteroclinic cycles of relevance to Eqs. (29-31). Structurally 
stable heteroclinic cycles (HETi) connecting two SR 3 solutions, out of phase by n (i.e., 
related by a translation Tu^^), have been well studied [34,36] in the context of the 1:2 
steady-state spatial resonance with 0(2) symmetry and in resonant triads as well [35]. They 
are generic solutions when af3 < that appear at small amplitude (|A|, <C 1). 

In addition, there are much more intricate heteroclinic cycles (HET2) involving connections 
among O, SR3, MSW K , and MM K (see [39]). Under appropriate conditions, these cycles form 
near the HETx cycles (in a neighboring region of parameter space but further from the 
origin). The associated dynamics is organized by a sequence of transitions among distinct 
heteroclinic cycles (one of which is structurally stable) and contains chaos of Shil'nikov type. 
The type of chaotic attractor one typically finds is shown in Fig. 6. 

When the resonance terms vanish (i.e., at a — (5 — 0) there are, under certain conditions 
[35], heteroclinic cycles connecting SRi, SR2, and SR3. Because of the resonant terms in 
Eqs. (29-31) these cycles are no longer present. Nonetheless, for large m + n when resonance 
effects are diminished, or at values of <p m and <p n where $ = ±7r/2 and the leading order 
resonant term vanishes (see Eqs. (32)), these "nearby" cycles may yet exert an influence on 
the dynamics. 

4-2 Bifurcations: numerical examples 

We provide here some numerical examples using coefficients derived from the Zhang- Vihals 
equations (41). Figure 7 shows the bifurcations that occur for m:n = 1:2, $ = 0, p = 1, 
v = 0.1, T = 21.5, and u = 200 (cgs units); for these parameters 7 = 0.0341 and we find 
a = -0.265, = 0.281, a = -6.136, b = 0.0266, c = -0.98, d = -2.076, and e = -2.234. 
We use m:n = 1:2 rather that mm = 3:2 (as in Section 3) mainly for pedagogical reasons: 
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Fig. 6. Chaotic attractor in Eqs. (50-53) obtained at (A,//) = (—0.0377,0.035) with 
the coefficients a = -0.19, (3 = 0.166, a = -0.0399, b = -0.436, c = 1.27, 
d = —1.43, and e = —0.106, derived from the Zhang- Vihals model (41) with 
m:n = 1:2, $ = tt/4, and (in cgs units) uj = 200, v = 0.01, T = 70. This at- 
tractor is in IA K and visits the neighborhood of two MM K fixed points: one with 
9 = and one with 9 = ir. The latter is near a Hopf bifurcation to MSW K . The 
boundary r% = should be identified at 9 = and 9 = tt. 

although the bifurcations seen with other forcing ratios such as 3:2 (and comparable physical 
parameters) are the same as those found in Fig. 7 they tend to be somewhat more bunched 
together thus more difficult to illuminate. The use of 1:2 forcing also has the advantage of 
the large 0(1/7) oscillations with <3> in b, c, and d (see Eqs. (32)); this leads to greater control 
over the normal form coefficients and hence to a wider range of dynamical possibilities (for 
example, this ^-dependence made it easier to find the interesting heteroclinic behavior shown 
in Fig. 6). 

If one traverses the clockwise circular path shown in Fig. 7 the resulting bifurcation diagram 
is that of Fig. 8. Initially, in the third quadrant of the (A,/i) plane, all initial conditions 
are attracted to the stable flat state O. As the polar angle (0, say) is decreased through 
7r this flat state becomes unstable to (a circle of) standing rolls of type SR3. At <fi ~ 2.64 
(point 1 in Fig. 8) these stable SR 3 states in turn undergo a supercritical symmetry-breaking 
bifurcation to (a torus of) /t-symmetric mixed modes MM K satisfying 9 = n. These MM K 
states are stable until a subcritical Hopf bifurcation at <fi ~ 1.64 (point 3) producing k- 
symmetric modulated standing waves MSW K . The MSW K are themselves destroyed in a 
heteroclinic bifurcation (point 2, marked by a small dash on the SR3 branch) joining O 
and SR 3 (see [39]) at <p ~ 1.67. The same global bifurcation marks the birth of structurally 
(and asymptotically) stable heteroclinic cycles (HETi) connecting points on the circle of 
SR 3 states to their u n translates". After the Hopf bifurcation renders MM K unstable the 
heteroclinic cycles HETi are the only attractors of the system (29-31). These HETi remain 
stable until ~ 1.16 (point 5; this is when fi = Xe/c and the two eigenvalues of Eq. (58), 
one positive and one negative, are of equal magnitude). The bifurcation that occurs here 
simultaneously produces stable MDW K and unstable MSW (the fact that the MSW branch 
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Fig. 7. Bifurcation sets for the case m:n = 1:2 and <3? = with (in cgs units) 
p = 1, v = 0.1, T = 21.5, and w = 200. Curves are labelled by the solutions 
generated (see text) and are dashed in the case of Hopf bifurcation. The clockwise 
circular path is used to generate the bifurcation diagram in Fig. 8 and is at a radius 
+ /x 2 = 0.0128. 
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Fig. 8. Bifurcation diagram for the clockwise circular path shown in Fig. 7 
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does not visibly originate here is due to numerical limitations). The stable MTW K branch 
describes modulations of ever decreasing amplitude as it approaches the branch of DW K 
states, themselves present since an 1Z symmetry-breaking bifurcation on MM K at ~ 1.63 
(point 4), and is destroyed in a Hopf bifurcation at ~ 0.95 (point 6). The DW K are then 
stable until a supercritical k symmetry-breaking bifurcation to DW at ~ —0.5 (point 8); 
this occurs just prior to their termination at ~ —0.55 (point 9) on the other MM K branch, 
i.e., the one with 9 = existing between = —tc/2 and a bifurcation from SR3 at ~ 0.29 
(point 7). The stable DW branch terminates, along with the unstable MSW branch, in a 
Hopf bifurcation on SRi at ~ —0.59 (point 10). The SRi rolls are stable between this 
Hopf bifurcation and a subcritical bifurcation to MM at ~ —0.79 (point 12). This MM 
branch also connects to (and stabilizes) the MM K branch at ~ —0.7 (point 11). Finally, 
the MM K states (with 9 = 0) are stable between this last bifurcation and = —n/ 2 where 
they vanish, along with SRi, restoring stability to the flat state. 

If, in contrast to the case of Fig. 7, one adjusts the forcing phases so that $ = n/2 then 
the (leading order) expressions for a and (3 of Eqs. (32) vanish. With the same physical 
parameters as Fig. 7 it turns out that the first nonvanishing contributions have the same 
algebraic sign, i.e., af3 > 0. Specifically, we find a = 0.0784, (3 = 0.0744, a = 0.549, 
b = -0.523, c = -0.0155, d = -3.266, and e = -2.227. The fact that a(3 > explains the 
dramatically different (and less rich) unfolding shown in Fig. 9. In this case there are no 



Fig. 9. Bifurcation sets for m:n = 1:2 with the same parameters as Fig. 7 except 
that $ = 7r/2. For this choice of phase a(3 > and there are no DW K , MDW K , DW, 
MSW K , MSW, or heteroclinic solutions. 

drifting solutions (DW K , MDW K , DW) or heteroclinic cycles (HETi, HET 2 ). Furthermore, 
there are no modulated standing waves MSW K or MSW (in was shown in [37] that af3 < 
is necessary for MSW K ). A second, somewhat surprising, feature of the $ = it/ 2 coefficients 



SRg 




MM„, 
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is the large change in the self-interaction coefficient a: it is now positive whereas at $ = it 
was —6.136. Such substantial ^-dependence is not an ingredient of Eqs. (32) and is a result 
of "1:2" linear resonances, i.e., harmonic modes satisfying k = 2ki or k = 2k 2 . The influence 
of these damped modes (see [10]) is especially strong here because k 3 ~ 1.86 is quite close to 
2&i ps 1.94 (the damping of the 2k\ modes is thus very small and they nearly satisfy the 1:2 
temporal resonance condition m = Q(2k\); in this case m = n/2). For comparison, note that 
the 2k% modes are not "almost" neutral and the corresponding self-interaction coefficient 
e shows no significant dependence on $, in agreement with Eqs. (32). Because a > the 
bifurcation of SRi and MM K from the flat state is now subcritical (observe how the MM 
bifurcations from SRi and MM K have switched to the left side of A = in Fig. 9) and the 
cubic truncation (29-31) is no longer well-behaved (some trajectories diverge). 



5 Conclusions 

The weakly nonlinear dynamics of a class of resonant triads occurring in the vicinity of the 
bicritical point of the two-frequency forced Faraday problem has been studied in detail. The 
approach taken for this resonant triad can be easily extended to a wide range of patterns 
produced by multi- frequency forcing (see [9]), and is based on very general symmetry consid- 
erations. In particular, we make use of the fact that the weakly damped problem possesses 
additional structure beyond the more obvious spatial and discrete temporal symmetries, 
i.e., it must satisfy the constraints of continuous time translation and time reversal (when 
these are viewed as parameter symmetries) and respect the nearby Hamiltonian structure 
of the undamped problem. The effect of these "broken" symmetries is straightforward to 
establish in the context of TW equations (17-18), and in the reduced SW equations (29-31) 
reveals itself in the structure of the normal form coefficients (32). Most notably, we found 
that the SW coefficients of the quadratic resonance terms (a and (3) could be expected to 
be of opposite sign, exhibit simple harmonic dependence on the phase $, and scale as the 
damping 7 to the power (m + n — 3)/2. Each of these predictions was supported in the direct 
numerical calculation of coefficients from the Zhang- Vinals equations (41). Moreover, the 
exponential decrease in the size of the resonance terms with m + n provides an explanation 
for the lack of experimental observations [3,8] of this triad with large values of m and n. 
Recall, however, that this simple prediction (as well as the other scaling laws) emerged after 
we expanded the SW normal form coefficients in powers of f m , f n , and 7 (see Eqs. (19-23)). 
Particularly in the case of 7 it is not clear that the assumption of analyticity is a good one 
- the addition of damping is known to act as a singular perturbation (see, e.g., [40,41]). For 
the quasipotential equations of Zhang and Vinals the numerical agreement (and hence the 
justification for a simple Taylor expansion) is quite convincing, but we speculate that a more 
rigorous treatment using the Navier-Stokes equations with realistic boundary conditions and 
finite fluid depth may lead to a departure, in some regimes, from the scalings given here. 
This departure will presumably be small in weakly damped systems of large depth where 
the Zhang- Vinals equations are expected to be valid. 

The resonant triad system (29-31) exhibits extremely rich dynamics, much more so in the 
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case where the resonant coefficients satisfy af3 < 0. That the underlying Hamiltonian struc- 
ture makes this the typical case for weak damping and, furthermore, that this situation 
depends on the relative phase of the two forcing terms, is a result we have emphasized 
throughout this paper. When a(5 < there exist several types of modulated standing waves, 
drifting waves (of both steady and modulated character), and heteroclinic cycles with associ- 
ated phenomena like the chaotic attractor of Fig. 6. It must be admitted, however, that none 
of these more complicated solutions has yet been reported in the experimental literature; 
other than simple roll states, only the mixed mode MM K has been established experimen- 
tally [8]. It is quite possible that the remaining solutions are simply unstable when all types 
of perturbations are considered (the system may prefer hexagons or squares, for example). 
It might also be that the use of finite containers plays a critical role in preventing many of 
these solutions from appearing (the drifting waves DW K , DW, and MTW K obviously cannot 
exist as such in nonperiodic bounded domains). It would be interesting to see if a more 
deliberate search, perhaps in a very large aspect ratio container, would reveal more features 
of the dynamics presented in Section 4. On the other hand, it would miss the point to focus 
too narrowly on the resonant triad treated here. This paper is also intended to serve as 
a detailed illustration of a very general method with wide applicability to parametrically 
forced problems, particularly those with forcing functions composed of multiple frequencies 
- for example, in [9] we considered two-frequency forced hexagons and two-mode superlat- 
tices [5]. The symmetry-based approach of this paper provides, for weakly damped systems, 
a means of determining scaling laws and other important qualitative features of the normal 
form coefficients without resorting to a full center manifold reduction from the governing 
equations (often an exceedingly difficult task). The results also suggest how one might be- 
gin to think about designing forcing functions that enhance, suppress, or otherwise control 
the patterns of interest (superlattices, for example). This very intriguing direction will be 
pursued in future work [42]. 
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A Appendix 



We give here the expressions for J ± , J 7 appearing in Eqs. (32). 

2g 2 r , v r \ n \ (■„,(', r„[>; 
Ji = —{Vn + ^ VnCx), where C x = — + H r§, (A.l) 



2^, ^ m , ,,. ' : [>„[>,. 

J<i = -jTT [Qm + — QnC 2 ) , where C 2 = — + — ~y + -7^-, (A.2) 



v 1 
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v m v r 


VnQr 
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QmVl 


! QnQr 


Qr 




ft. A? 


~2 ' 



J 3 = a r + b r - Ci(aj + bj) + -^((^j - (qi)j - (q 2 )t), (A. 3) 



A 



■u. 



J 4 = S r + c r + d r - dfa + Ci + di) + -^((rs); - (ri)i - (q 3 )i - (q 4 )i), (A.4) 

Aj 

J 5 = e r + f r - Ci(ej + fi) - Y"((q5)i + (qe)i) + -=-({ T s)i - (r^i), (A.5) 

A« A 4 * 

J 6 = g r + h r - C 2 ( gi + \) - ^(( Sl ), + (s 2 ),) + ^((u 2 ), - (ui)i), (A.6) 

A*i Aj 

J 7 = l r + Pr - C7 2 (li + Pi) + ^ ((U3)i - (s 3 )i - (s 4 )i) • (A.7) 
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